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Abstract.  The  linear  stability  analysis  of  the  / = 1 diocotron  perturbations  in  a single  charged 
plasma  confined  in  a cylindrical  Penning  trap  is  critically  revisited.  Particular  attention  is  devoted  to 
the  instability  due  to  the  presence  of  stationary  points  in  the  radial  profile  of  the  azimuthal  rotation 
frequency.  The  asymptotic  analysis  of  Smith  and  Rosenbluth  [1]  for  the  case  of  a single-bounded 
plasma  column  (algebraic  instability  proportional  to  f1/2)  is  extended  to  the  case  of  a cylindrical 
Penning  trap  with  an  additional  coaxial  inner  conductor,  and  it  is  shown  that  the  algebraic  instability 
found  in  the  case  of  a single-bounded  plasma  column  becomes  exponential  at  longer  times.  The 
relevant  linear  growth  rate  is  computed  by  a suitable  inverse  Laplace  transform  (contour  integral  in 
the  complex  plane).  The  analytical  results  are  compared  with  the  numerical  solution  of  the  linearized 
two-dimensional  drift  Poisson  equations. 


INTRODUCTION 

The  linear  stability  analysis  of  the  Z = 1 diocotron  perturbations  in  a low  density  single 
charged  plasma,  radially  bounded  by  two  cylindrical  conductors  held  at  fixed  potential, 
is  critically  reviewed.  Using  a model  with  a radial  step  density  profile,  Levy  [2]  showed 
that  the  plasma  is  neutrally  stable  when  it  is  in  contact  with  one  or  both  conductors,  or  if 
the  charge  on  the  inner  conductor  is  large  enough.  If  the  central  conductor  is  absent,  the 
l = 1 diocotron  mode  is  neutrally  stable,  while  lower  l > 2 modes  may  be  unstable.  The 
effect  of  a central  conductor  on  the  stability  of  an  hollow  plasma  column  has  been  also 
studied  experimentally  [3]. 

Particular  attention  is  devoted  here  to  the  instability  due  to  the  presence  of  one  or 
several  stationary  points  in  the  radial  profile  of  the  azimuthal  rotation  frequency  of  the 
plasma.  The  asymptotic  analysis  of  Smith  and  Rosenbluth  [1]  for  the  case  of  a single- 
bounded  plasma  column  (instability  proportional  to  r1/2)  is  extended  to  include  algebraic 
instabilities  growing  as  f“,  with  1/2  < a < 1.  The  asymptotic  analysis  is  generalized  to 
the  case  of  a trap  with  a coaxial  cylindrical  inner  conductor,  and  it  is  shown  that  the 
algebraic  instability  found  in  the  previous  case  becomes  exponential  at  longer  times:  the 
relevant  linear  growth  rate  are  computed.  Finite  length  and  finite  Larmor  radius  effects 
are  neglected. 
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FIGURE  1.  Scheme  of  a cylindrical  Penning  (rap.  Loft:  without  central  conductor;  Right:  with  central 
conductor. 


BASIC  EQUATIONS 


In  the  model,  a one  component  plasma  is  assumed  to  be  contained  within  two  infinitely 
long  cylindrical  conducting  walls,  of  radii  a and  /?,  respectively.  The  external  conductor 
is  grounded  (see  Fig.  1 ).  The  case  where  the  central  conductor  is  absent  is  treated  simply 
by  setting  a = 0.  The  system  is  immersed  in  a static  and  uniform  magnetic  field  B = Bz. , 
directed  along  the  axis  of  the  trap.  In  the  model  considered  here,  the  evolution  of  the 
system  is  described  by  the  two-dimensional  drift  Poisson  equations  [4],  written  in  polar 
coordinates  (r,0), 

^ + (V<px  V/i)  z = 0,  V2q > = /!.  (1) 

Adimensional  quantities  are  used.  The  density,  n,  is  normalized  over  an  arbitrary  ref- 
erence density,  n;  the  lengths  over  the  radius  of  the  outer  conductor,  r = R (an  explicit 
notation  for  R is  kept  in  the  following);  the  frequencies  over  cb  = 4ne2n/m(i>c,  where 
(0C  = -eB/mc  (e  and  m being  the  charge  and  the  mass  of  the  particles,  respectively, 
and  c the  velocity  of  light  in  vacuo);  the  potential,  (p,  over  (p  = 4nenr2;  and  the  electric 
charge  per  unit  length  on  the  central  conductor,  Q,  over  Q = Kefir2. 

Linearizing  Eqs.  (1)  for  perturbations  6cp  = <P/(r,/)e',e  and  8 n = »/(r,/)e'/e  with  a 
given  azimuthal  number  /,  yields  the  following  second-order  differential  equation  for 
the  potential  amplitude  <)>/  (see,  e.g.,  Ref.  [4]): 


ra  , i 

i a d i2' 

— + t/O)f(0 

r dr  'dr  r2 

<W  -jn'o(r)h  = 0; 


1 f 

<%('')  = 3 / >H)(r)rdr  + 
I J a 


Q 

2r2 


(2) 


where  «o(r)  and  (0£-(r)  are  the  unperturbed  density  and  the  unperturbed  azimuthal 
frequency  of  the  plasma,  respectively,  and  a prime  denotes  the  radial  derivative.  Eq.  (2) 
has  to  be  supplemented  with  the  boundary  conditions  |<)>/(0)|  < +<»,  <)>/(/?)  = 0 in  the  case 
without  central  conductor,  and  <p/(«)  = 0,  $/(/?)  = 0 in  the  case  with  central  conductor, 
respectively. 
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In  the  following,  the  analysis  is  restricted  to  the  case  1=1.  The  Laplace  transforms 
for  the  perturbed  potential  and  density  can  be  written  in  this  case  as 


fK  hi 

<Mr)  = rlp+ '“eO)] Jr 

in'o(r) 


R h(x)-h(p) 


nP(r)  = 


p + mK(r) 


\p  + mE  (a-)]2  dx’ 

R h(x)  — h(p) 


X3[p+i(0E(x)]2 


dx, 


where 


h(r)  = f x2ni(x,0)dx;  h(p)  = f 

•la  Ja  X 


h(x)dx 


[p  + /a>E(.x)]2 


f 

J a 


dx 


x3  [p  + i'0)e(x)]: 


(3a) 

(3b) 

■■  (4) 


The  inversion  of  the  Laplace  transforms,  e.g.,  $](/',/)  = / fe<| )p(r)e,n,  involves  an 

Br 

integration  in  the  complex  plane  along  a Bromwich  contour,  Br,  which  goes  to  the  right 
of  all  the  transform’s  singular  points  and  branch  cuts. 


TRAP  WITHOUT  CENTRAL  CONDUCTOR 


The  inverse  Laplace  transformation  can  be  performed  in  a closed  form  if  a = 0,  recov- 
ering the  result  of  Ref.  [1]: 


<E>]  (>',?)  = r J [/i(x)/x3]  [1  + /(0£(r)/-/CO£(x)t]  e ,0)£W'  dx; 
nt(r,t)  = n , (r, 0)e-/(O£W'  - itn'0{r)  j\h(x)/x3]e^i(,iE(x)l dx. 


(5a) 

(5b) 


It  is  readily  seen  that  an  unstable  perturbation  can  not  grow  faster  than  t,  as  the  integrals 
in  Eqs.  (5)  are  decreasing  functions  of  time.  An  asymptotically  non-decaying  solution 
can  originate  from  the  density  jumps,  from  the  two  ends  of  the  interval  of  integration  in 
Eqs.  (5),  or  from  the  stationary  points  inside  the  same  interval. 

If  n(r)  is  non-monotonic,  (0£(r)  may  also  be  non-monotonic  with  its  extrema  being 
the  points  of  stationary  phase  for  the  integrand  functions  in  Eqs.  (5).  If  r = ro  is  a 
nondegenerate  stationary  point  for  o%  (tOg(ro)  = 0,  co'^(ro)  ^ 0),  then 


<t>i 


i\Z2nrh(ro) 

V/I®£(ro)ko 


[co£(r)  - 0D£(ro)]  H(r0  -r)\ft  exp 


in 


-ito£(ro)t- jsignjfflg) 


(6) 


where  H(x)  is  the  Heaviside’s  step  function.  The  perturbation  grows  in  this  case  propor- 
tionally to  sft  [1]  (SR-instability). 

In  the  presence  of  an  inflection  point,  r*,  of  the  (0 E profile,  different  situations  may 
occur,  according  to  the  value  of  (^(r*)/®^  (r*).  If  a/E(rt)/(0^ (r*)  < 0,  0)g(r)  has  two 

extrema  at  r±  = r*  ± ^ -2p)'E / (pg . At  large  times,  t » |r|_1,  with  T = 
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FIGURE  2.  Left:  Smooth  h«  profile  (solid  line)  and  corresponding  (0£  profile  (dashed  line),  which  give 
rise  to  exponentially  damping  initial  perturbations.  Right:  Amplitude  of  the  perturbed  potential,  |4>i  |(r,t), 
vs.  r and  f. 


these  extrema  contribute  separately  to  <J>|,  and  the  initial  perturbation  grows  as  \ft.  If 
C0/.  (r,)  = 0,  these  points  merge  at  r = r,,  and  the  perturbation  grows  proportionally  to 


,2/3 ; 


<t>l  ~ 


t'2|/-1r(l/3)r/t(r,) 

3'/6|a)»'|'/3,-3 


[cot- (r)  - CO t-(rt)] H(r,  - r)t2/3e 


(7) 


Finally,  if  (o'£(r*)/o)£  (r»)  > 0,  the  extremum  points  move  into  the  complex  plane,  while 
4>i  grows  algebraically  according  to  Eq.  (7)  in  the  initial  time  evolution,  t < |T|  '.and 
then  decays  exponentially  at  late  times  (this  situation  is  illustrated  in  Fig.  2): 


4*i  ~ 


iV2irh(rt) 
(2«),/4  rj 


[cot(c')  - C0t(r, )]  H (r»  -r)\ft  exp[-iffl/T  - Tf] . 


(8) 


The  frequency  of  the  “saturated”  mode  is  c Df(/?)  (Levy’s  mode).  In  general,  if  the 
a >t  profile  presents  a stationary  point,  /*,  of  order  m (c o\P(rt)  = 0 for  j < m , and 
co^(r„)  7^  0),  the  following  asymptotic  behavior  of  the  perturbed  potential  is  obtained: 


<t>l  ~ 


2 

m 


m\ 

4”V.) 


r* 
’ * 


[co£(r)  - CO/.-  (r*)]  H(rt 


r) 


'(-I)"'-1  - r 
2 


sign^'V*))  sin(^)}- 


(9) 


For  m = 2,  this  formula  reduces  to  Eq.  (6),  and  for  m = 3 to  Eq.  (7). 
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Finally,  a linear  growth  of  the  amplitude  of  an  initial  perturbation  is  found  if  an 
interval  exists  where  ©£(r)  is  constant,  and  if  this  interval  does  not  contain  the  whole 
plasma  column.  This  situation  can  be  realized  by  means  of  a two-column  plasma,  where 
the  densities  of  the  internal  (0  < r < bo)  and  external  (a \ < r < b\)  columns  are  related 
by  tio(bo)  b^  = no(a\ ) «2.  In  this  case: 

<|)|  w itr[mE(r) -a)E(a\)]H(a\ - r)  e~'“c(fll )'  f [h(x)/xi]dx.  (10) 

Ja\ 


TRAP  WITH  CENTRAL  CONDUCTOR 


The  SR-instability  comes  from  the  neighborhood  of  a nondegenerate  stationary  point 
ro  of  ©£.  This  instability  can  be  interpreted  as  the  contribution  of  the  branching  point 
p =z  — /co/j(r0)  in  the  complex  p-plane.  When  a central  conductor  is  present  in  the  trap, 
it  results 


h(p)  ~ 


Jf(p  + /(0E(ro))/>-3 
J(p  + i(0E(ro))/rl  + fi{p)  ,0 


(11) 


ne-/jtsign[(O^(r0)]/4  j 

i/j2ai|(ro)|o3/2  ’ 2fl2[p  + iA©£][p  + !©£(a)]  ’ 


where  A(0£  = no(a)/2  and  ©e(u)  = Q/2a 2,  respectively.  The  function  J(a)  is  two- 
valued since  it  contains  the  rational  power  3/2  of  a = p+ia>E(ro).  The  “physical  sheet” 
of  the  Riemann  surface,  arg[a]  € ( — Tt, 7t),  corresponds  to  the  complex  a-plane  with  the 
branch  cut  a € (-°°,0].  The  effect  of  the  central  conductor  turns  out  to  be  negligible  at 

small  times,  t « T « Ifflgl1/3©/4^3^-4/3,  while  at  larger  times,  t » T,  the  algebraic 
SR-instability  disappears.  The  inverse  Laplace  transform  can  be  performed  (asymp- 
totically for  t ->  <»),  by  means  of  a suitable  deformation  of  the  inversion  (Bromwich) 
contour  inside  the  “physical  sheet”.  In  particular,  an  exponential  instability  is  found  if 
[a>£(ro)  — A©£][co£(ro)  - ©£(a)]  > 0,  with  a growth  rate  (in  the  limit  a — > 0) 


ol/2_2/3  4/3 

^ 92/3  2 TrV^Tl/3  E ^ ~ [C°£(ro)  “ ■ 

z / r0  |©£|  i 


(12) 


Typical  results  of  the  analysis  and  the  numerical  simulations  of  Eqs.  (2)  in  the  case  with 
central  conductor  are  shown  in  Figs.  3-4. 


CONCLUSIONS 

The  existence  of  / = 1 diocotron  instabilities  in  a charged  plasma  confined  in  a cylin- 
drical Penning  trap,  growing  with  time  faster  than  the  SR-instability  (°c  ta , with  1 /2  < 
a < 1 ),  has  been  pointed  out.  In  addition,  it  has  been  shown  that  the  presence  of  an  inner 
conductor  (even  very  thin  and  uncharged)  can  transform  the  algebraic  instability  into  an 
exponential  one  at  late  times.  A criterion  for  the  occurrence  of  the  exponential  instability 


321 


*Jf).  2ajr) 


FIGURE  3.  Left:  «o(r)  (solid  line)  and  cd£(/-)  (dashed  line)  profiles,  which  determine  a SR  algebraic 
instability  [1]:  no  = [l  + {r/rp)2/ A]  [l  - (// r;,)2]‘  for  r < rp,  n«  = 0 for  r > r,„  with  A = 0.25,  rp  = 0.6. 
Center:  Amplitude  of  the  perturbed  potential.  |$i|(r,r),  vs.  r and  t.  Right:  Amplitude  of  the  perturbed 
potential,  |4>i  ((r, r),  vs.  r and  /,  for  the  density  profile  truncated  at  the  radius  a = 0.05  of  the  central 
conductor.  Q corresponds  to  the  total  charge  of  the  particles  lying  within  r = 0.05  in  the  previous  case. 


FIGURE  4.  Left:  Generalized  parabolic  density  profile  (solid  line)  used  in  the  computations:  ;in  = 
[I  — (r  — r,)2/r2]2  for  |r  — ;y|  < rp  and  /»,.  = 0 otherwise  (with  r,  = 0.5  and  rp  = 0.1)  and  co/;  profiles 
(dotted  lines),  plotted  for  a = 0. 1 and  different  values  of  Q.  Right:  Example  of  frequency  spectrum,  for 
£>  = 0 and  a = 0.2,  0.1, 0.05,  0.025,  and  0.001 . The  continuum  spectrum  lies  on  the  Re[co]  axis  from  0 to 
to£-(ro).  For  a given  a,  there  arc  two  complex  conjugate  discrete  frequencies,  which  converge  to  ro//(»o)  as 
a 0.  Stable  solutions  (Im[o)]  < 0)  do  not  belong  to  the  “physical”  sheet  of  the  Ricmann  surface. 


and  the  growth  rate  have  been  computed  analytically  in  the  limit  of  thin  inner  conduc- 
tor. This  work  has  been  supported  by  the  Italian  Ministry  of  Education  and  Scientific 
Research. 
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